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ABSTRACT 


The method of multiple scales Is used to derive a nonlinear 
Schrddinger equation for the temporal and spatial modulation of 
the amplitudes and the phases of waves propagating in a hard- 
walled circular duct. This equation 1s used to show that mono- 
chromatic waves are stable and to determine the amplitude dependence 
of the cut-off frequencies. 
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Introduction 


Nonlinear wave propagation In hard-walled ducts was Investigated 
by Coppens^, Pestorlus and Blackstock*, Maslen and Moore*, Burns'*, 

Keller and Mlllman*, and Keller*. Coppens^ and Pestorlus and Blackstock* 
determined the viscous and thermal dissipative effects on the nonlinear 
propagation of plane waves. Maslen and Moore* used the method of 
strained parameters (e.g., Sec. 3.1 of Ref. 7) to analyze strong trans- 
verse waves In a circular cylinder. Burns'* obtained a straightforward 
expansion limited to small axial distances because It contains secular 
terms. Keller and Mlllman* determined the nonlinear wavenumber shift 
of the symmetric dispersive modes by using the method of strained 
parameters. Keller* determined the amplitude dependence of the cut-off 
frequencies for the case studied by Keller and Ml liman*. 

The nonlinear effects of acoustic materials were investigated 
experimentally by Zorumski and Parrott® and Kurze and Allen* and 
analytically by Kurze and Allen*, Ingard*®, Isakovich**, and Nayfeh and 
Tsai**’**. Nayfeh and Tsai***’*® analyzed the combined nonlinear effects 
of the gas and the acoustic material. 

The purpose of the present paper Is to analyze the nonlinear pro- 
pagation of a wave packet by deriving a nonlinear Schrbdinger equation 
governing the temporal and spatial modulations of the amplitudes and 
the phases. 
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I. PROBLEM FORMULATION 


We consider finite-amplitude waves propagating In a hard-walled 
cylindrical duct. Although the analysis Is valid for any non-rectai^ular 
duct cross section, wa treat the case of a circular cross section so that 
we would be able to give an explicit solution. The fluid Is assumed to 
be Inviscid, irrotational, and Initially quiescent with a uniform pres- 
sure pt and a uniform density pt so that its subsequent motion can be 
represented by a potential function. 

We Introduce a cylindrical coordinate system (r, 6, x) whose x axis 
coincides with the duct axis. Dimensionless variables are Introduced 
by using the radius of the duct R*, the ambient speed of sound c|, and 
the ambient density p* as reference quantities. Thus, we let ? ® r*/R*. 

V e v*/ctt p * p*/p?, p = pf/ptet*. and t » t*ct/R*, where starred and 
unstarred quantities denote dimensional and dimensionless quantities, 
respectively, r is the position vector, v Is the velocity vector, p Is 
the gas density, p is the gas pressure, and t Is the time. In terms of 
these dimensionless quantities, the equations describing the conservation 
of mass and momentum are 

If + V • (pv) » 0 (1) 

p(|f + V »7v) « - 7p (2) 

The pressure Is related to the density by the Isentropic relationship 

P*/pt “ (p*/Pt)^ 

or In dimensionless quantities by 

YP - P^ (3) 
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where y is the gas specific heat ratio. Since the duct walls are 
assumed to be rigid, the appropriate boundary condition is the vanishing 
of the normal velocity at the duct walls; that is, 

V = 0 at r = 1 (4) 

where v is the radial component of velocity. 

Since the flow is assumed to be inviscid and irrotational , the 
velocity v is derivable from a potential function (j){r,t) according to 

V = V4> (5) 

Substituting for p and v from Eqs. 3 and 5 into Eq. 2, using the 
irrotational ity of the gas, and integrating the resulting equation, we 
obtain 

pT . 1 + {1 . y) (6) 

Eliminating p from Eqs. 3 and 6 gives 

/ 1 \y/(y-i) 

YP = (l + (1 - Y)[<l>t + ^V4.)23j (7) 

Differentiating Eq. 6 with respect to t, eliminating p by using Eqs. 1, 

5, and 6, and arranging, we obtain 

- V*(|> = (1 - Y)C<J>t + 

- 1^ (V(j.)2 - i74»-V(V4>)* (8) 

In terms of the potential function 4», the boundary condition 4 becomes 

* 0 at r * 1 (9) 

To determine an approximate solution to Eq. 8 subject to the boundary 
condition 9, we use the method of multiple scales (e.g., chapter 6 of 
Ref. 7 ) and let 
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♦(r.e.x.t) - I cV(r. 0 ,Xo,X,,X 2 . ToJi.Tj) + O(e-) (10) 

n«i " 

where e is a small but finite dimensionless parameter characterizing the 
amplitude of the wave and 

Xn » e"x, » e^t (11) 

Here, X# is a short scale characterizing the wavelength, Xi and X 2 are 
long scales characterizing the amplitude and phase modulations with 
axial distance. To is a short scale characterizing the frequency of the 
wave, and Tj and T 2 are long scales characterizing the temporal ampli- 
tude and phase modulations. Using Eq. 11 and the chain rule, we express 
the temporal and axial derivatives as 


9 _ 9 4 . ^ 9 4 . 

( 12 a) 

3-3 J.C.3 j.^2 3 . 

at- 5T7* ' sTT ' aTT • ' • 

( 12 b) 


Substituting Eqs. 10-12 into Eqs. 8 and 9 and equating coefficients of like 
powers of e, we obtain 
Order e 



(13«) 

„2 - 3 * .13 .1 3* , 3^ 

03b) 

34»i/3r » 0 at r » 1 

04) 


Order 

^ (♦.> • Ir^lTr * 2 Ife: * <1 - y) |f^ v?*. - ds) 

34>2/3r ■ 0 at r ® 1 (16) 
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Order e* 




- 2 



8^4»i 

3Xo3X2 


* ^ Ift 

+ (1 - y) ^o4'i'*’ 2(1 ■ y) (^o4>j 

-2 -11^) - 2 (Vo$i •Vo'J’z) - ■^ Vo<|)i *7o (Vo4>i)* (17) 

3(|)3/3r = 0 at r = 1 (18) 

II. SOLUTION 

We take the solution of Eqs. 13 that is bounded at the axis in the 
form of a traveling wave packet centered at the frequency w and the wave- 
number k; that is, we let 

«j>i= A(Xi,X 2 ,Ti,T 2 )Jjjj(<r) exp (i;) + cc (19a) 

? = kXfl - ujTo + m9 (19b) 

where J is Bessel's function of order m, cc stands for the complex 
m 

conjugate of the proceeding terms, and u and k satisfy the dispersion 
relationship 


w* - k* * K* 

(20) 

Substituting Eqs. 19 into Eq. 14, we have 

j'(.) . 0 

(21) 
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In what follows, we exclude the non-dlspersive case corresponding to 
K > 0 (I.e., plane waves). Note that the function A Is still undeter> 
mined at this level of approximation; it 1s determined by Invoking the 
so-called solvability condition In the second-and third-order problems. 
Substituting for (|)i from Eq. 19a into Eq. 15, we obtain 

-/ ( 4 ) 2 ) = 21(o) + k |^)Jj^(Kr)exp(ic) 

+ [21a)K*J^^(icr) - 4iu)<m f 

- 1a)(YtD* - + k^) J|^(<r)]A^ exp(2ic) + cc (22) 


Since the homogeruous second-order problem consisting of Eqs. 16 and 
P2 is the same as the first-order problem and since the latter has a 
non- trivial solution, the inhomogeneous second-order problem has a 
solution if, and only if, a solvability condition is satisfied; this 
condition yields 


With this solvability condition, the solution of the second-order 
problem is 

^ [riJ^(icr) + ^ r3J2j^(2Kr)]A^ exp(2ic) + cc 


where 


(2fia) 


Fi = - iu[j m(y + l)o)^ic”® + 1], F 2 ®^i(y+ l)oj^< *, 






(24b) 
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Substituting for and i ))2 from Eqs. 19 and 24 Into Eq. 17, we 

have 

JT Us) - [21(w |^+ k 1^) • |r|+ 

- A*AF(r)exp(1c) + cc + harmonics other than exp (tic) (25) 

where F(r) is given In Appendix A. Since the homogeneous third-order 
problem consisting of Eqs. 18 and 25 has a non-trivlal solution, the 
corresponding inhomogeneous problem has a solution If, and only If, 
a solvability condition Is satisfied. To detemlne this solvability 
condition, we seek a particular solution of the form 

♦s * tj;(r,Xi,X2,Ti,T2)exp(1c) (26) 

Substituting this solution Into Eqs. 18 and 25 and equating the co- 
efficients of exn(lc) on both sides, we obtain 

* Irf * 

3tj)/3r ■ 0 at r * 1 (28) 

Multiplying Eq. 27 by rJjjj(tcr), Integrating by parts from r » 0 to r * 1, 
and using Eq. 28. we obtain the solvability condition in the form 

It 7 ^ ‘ Irf Ixf “ ^ 

where 

A « [} rF(r)J (Kr)dr]/[} rJ^(Kr)dr] (30) 

III 0 III 

Eliminating 3^A/dXi from Eq. 29 by using Eq. 23 gives 

21(u) |t 7 + k 11^) + (|r - 1) Ify *A A^^ (31) 
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To simplify Eq. 31, we differentiate the dispersion relationship 20 
with respect to u and obtain 

kk* « u) (32) 

where k" = dk/du the Inverse of the group velocity. Differentiating 
Eq. 32 with respect to u gives 

kk" - 1 . k'* - 1 - wVk* (33) 

Using Eqs. 32 and 33, letting » e”t and Xj ■ e*x, and arranging, 
we rewrite Eq. 31 as 

(34) 

Changing tho Independent variables from x and t to 

E ■ t - k*x, n “ X (35) 

we express Eq. 34 In the form 

(36) 

which Is a nonlinear Schrodinger equation. Letting A = j a exp (1$) with 
real a and 6 In Eq. 36 and separating real and Imaginary parts, we obtain 

3^ • '< Lh 35 * 7* 3^J ° 


li 

3n 


+ 





(38) 


The Case of Monochromatic Waves 

For monochranatic waves, 3a/3C *36/3C • 0, and Eqs. 37 and 38 can 
be Integrated to give 


9 


(39) 


• • •* «nd e ■ - g- e* ^ ajn ♦ Bo 

whert ao and Bo art constants. Substituting for a and 6 from Eqs. 39 
Into Eqs. 19 and 24, we rewrite Eq. 10 In terms of the original 
variables as 


^ ■ eaJ|jj(icr) cos (ibc - wt me ♦Bo) 

♦ ‘ [ r,J>(Kr) ♦ rjrJ„ (icr) (ter) * r,J,„ (2n-)j 

/ 

COS [ 2(tb( - <i)t ♦ roe ♦ 3p) ] + 0(c*) (40) 

where 

l( • k - J c*Ak'*aj (41) 

This result can be obtained as a special case from the solution of 
Nayfeh and Tsai** by letting the resistivity of the acoustic material to 
be Infinite. Moreover, our solution reduces when m « 0 (I.e., symmetric 
modes) to that of Keller and Ml'ilman, which they obtained using the 
method of strained parameters. 
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Eqi ‘.Ion 41 shows that the nonlinearity of the gas results In a 
shift In the waven»anber from k to Ic ■ k - je*Ak->a' The ruaaerlcal 
results of Nayfeh and Tsai*' show that the nonlinearity shifts the 
wavenumber to lower values resulting In higher phase speeds. The 
wavenumber shift Increases with Increasing frequency and decreases 
with Increasing azimuthal or radial mode manber. 

Equations 37 and 38 can be used to analyze the stability of the 
aforementioned monochromatic solution. To do this, we let 

a > ao ‘t- ai , 8 ■ - ♦ 8o S: (42) 

where ai and 8t are small compared with the preceeding terms. Sub- 
stituting Eq. 42 Into Eqs. 37 and 38 and neglecting the nonlinear terms 
In ai and 6i, we obtain 



(43) 


38^ + i It'*- JL- 


3n 


So 




A * * 
Ic *•** 


(44) 


Since Eqs, 43 and 44 linear, we seek their solution In the form 

ai « Ij exp£l(kn - iSE)3, 8i “6i expf1(k n - wE)] (45) 

where Si and 8i are constants. Substituting this solution Into Eqs. 43 
and 44 and eliminating ii and Bi , we obtain 

it* • j k"*o5* (a* - j e*Aa|/kk") (46) 

which shows that. If A/k"<0, k Is always real for all values of u 
so that the monochronmtic waves given by Eqs. 40 and 41 are neutrally 
stable. On the other hand. If A/k**>0, k* Is negative for all 
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w < ca» (A/2kk")*/*; consequently, disturbances grow exponentially 
with C and monochromatic waves are unstable. Since the numerical 
results of Nayfeh and Tsai** show that A Is positive, monochromatic 
waves are sUble only if k" < 0. However. Eqs. 20 and 33 show that 
k“ » -icVk* <0. Therefore, monochromatic waves are stable. 

Solution Near Cut-Off Frequencies 

Although Eqs. 40 and 41 are valid for a wide range of frequencl3s, 
they break down as k-*^0 (i.e., near the linear cut-off frequencies) 
because the wavenun4>er shift approaches infinity. However, the basic 
equations and 29 for the modulation of the amplitude and the phase 
with axifl di<,tance and time are valid for all frequencies. In this 
section, we specialize these equations to frequencies near the cut-off 
values. To do this, we use Eq. 23 to express 3A/3T| as -(k/u) dA/3X}, 


substitute this result into Eq. 29, and obtain 

^ ^ A**" (47) 

Letting Ti « e*t and • c^x in Eq. 47, we rewrite it as 

Jiu ♦ 21k ^ * 0 - y) 0 ■ t’A A’f (4Sa) 

ll * lx ■ 7 If ■ ■ A'S (48b) 

For monochromatic waves, dA/at > 0 and Eq. 48a becomes 

«k ^ ^ (1 - ^) 0 - a'A ATT (49) 

which is valid for all frequencies away from zero. 
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( 50 ) 


Equat'on 49 has solutions of the form 

A ■ |a exp (ie) 

where a Is constant and 

( - k ♦ Ck’ - JcMl - ^) (1 - ^)-> (51) 

Away from the cut-off frequencies, k is away from zero and the radical 
in Eq. 51 can be expanded for small c yielding 

II = - J k-> (52) 

in agreement with the monochromatic solution obtained above. On the 
other hand, when k-K) (i.e., near the cut-off frequencies), Eq. 51 tends 
to 

11= - k+ (k* - }e'Aa>)‘'* (53) 

Substituting for A from Eq. 50 into Eqs. 19 and 24, using Eq. 53, 
and letting To = t and = e'^x, we obtain Eq. 40; however, ic of Eq. 41 
is modified to 

k = (k" (54) 

Therefore, the cut-off frequencies are solutions of 

k^ - 1 = 0 (55) 

Since k* = to* - according to Eq. 20, the cut-off frequencies are 

u) = ic + ^ e*a^A (kc) + . . . (56) 

where A (ic) stands for the value of A when oj = ic. These cut-off 
frequencies reduce when m = 0 to those obtained by Keller*. 
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Using Eq. 48a, one can carry out a stability analysis and show 
that this modified solution Is stable. 
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Appendix A 

F(r) = k\M j;*(Kr) [^y-D <*>2 + - 3k^ + 2k*J 

- K Jj (Kr) (icr) + 2iu + V 2 <r)^ 

^ '^tn+i 

[- ■ 2riKJm(Kr) J r3J^^(2Kr) - 2r3K)^,^^ 

- [ (y- 1V* + ^ + 2k"] [|^+ + (y-1)u 3" + 2iu)r,] 

+ 2i« J^(Kf) [r^YJ^CKr) J^+j(<r) + T3 J,^(2icr)]]y 
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